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ELEMENTS  AND  PARAMETERS  OF  THE  OSCULATING  ORBIT 

AND  THEIR  DERIVATIVES- - 

By  Wilbur  F.  Dobson,  Vearl  N.  Huff, 
and  Arthur  V .  Zimmerman 


SUMMARY 

The  analysis  determines  the  time  derivatives  of  the  conic-section 
orbital  elements  in  a  disturbed  orbit  by  perturbation  methods.  Integra¬ 
tion  of  any  of  the  several  resulting  systems  of  six  first-order  linear 
differential  equations  by  numerical  methods  can  be  a  useful  tool  for  the 
solution  of  problems  in  orbital  mechanics. 

Equations  for  the  two-body  orbit  are  also  summarized  in  a  convenient 
manner . 


INTRODUCTION 

The  recent  emphasis  on  space-flight  trajectory  calculations  and  the 
use  of  electronic  computing  machinery  have  combined  to  increase  the  in¬ 
terest  in  the  perturbation  methods  for  studying  problems  in  celestial 
mechanics.  This  is  especially  true  of  the  methods  that  leave  the  dis¬ 
turbing  function  undeveloped  and  require  numerical  integration  to  com¬ 
plete  the  solution.  Previously,  the  amount  of  numerical  work  required 
rendered  precision  calculations  by  these  latter  methods  impractical.  It 
is  the  purpose  of  this  report  to  examine  and  extend  some  of  the  previous 
work  in  perturbation  theory  to  secure  fonas  that  may  be  better  adapted  to 
numerical  integration,  at  least  for  specific  problems . 

The  basic  work  in  developing  expressions  for  the  derivatives  of  the 
orbital  elements  must  be  credited  to  various  notable  contributors  in 
dynamical  astronomy.  Perturbation  theory  was  begun  by  Euler'  in  1748. 
However,  the  first  complete  developnent  was  presented  by  Lagrange  in 
1782. 


The  perturbation  method  summarized  herein  is  formulated  in  terms  of 
Lagrangian  brackets.  Numerous  methods  for  evaluating  the  brackets  have 
been  published.  The  indirect  method  of  evaluating  the  brackets  used 
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herein  is  attributed  to  Whittaker^  as  reported  by  Smart  (ref.  l) .  The 
characteristic  of  the  indirect  methods  is  that  th^  work  begins  with  the 
derivation  of  a  general  expression  for  a  Lagrangian  bracket^  from  which 
all  brackets  are  easily  evaluated. 

For  convenience,  the  present  report  gives  alternative  forms  of  the 
perturbation  derivatives;  and,  by  using  the  results  presented  in  the 
tables,  many  others  are  obtainable.  The  extension  to  the  case  of  circu¬ 
lar  orbits  has  been  included.  A  collection  of  useful  two-body  equations 
is  also  given  in  table  I  without  derivation. 

The  proced\are  indicated  for  the  reduction  of  the  three  second-order 
linear  differential  equations  of  motion  in  rectangular  coordinates  to 
six  first-order  linear  differential  equations  in  orbital  elements  follows 
the  pattern,  but  is  revised  from  that  given  in  Moulton  (ref.  2).  The 
analysis  has  been  further  generalized  by  avoiding  the  requirement  that 
the  pertiarbation  fimction  be  a  potential  function.  This  extension  shows 
that  the  resiolts  are  valid  for  thrust  and  drag,  which  are  not  potential- 
type  functions.  Another  revision  concerns  the  determination  of  the  dis¬ 
turbing  functions  in  terms  of  the  elements  that  define  the  size,  shape, 
and  position  in  the  orbit.  This  procedure  ^iven  herein  is  believed  to 
be  more  direct  than  that  given  elsewhere  (refs.  1  and  2). 


ANALYSIS 


Equations  of  Motion 


Consider  the  motion  of  an  object  subject  to  an  inverse-square  cen¬ 
tral  gravitational  acceleration  directed  toward  the  origin,  and  also  sub¬ 
ject  to  smaller  disturbing  accelerations  that  can  be  expressed  as  func¬ 
tions  of  the  variables  and  constants  of  the  problem.  Let  OX,  OY,  OZ  in 
figure  1  be  the  coordinate  axes  in  a  noninertial  Cartesian  system  having 
its  origin  located  at  the  center  of  the  mass  M^.  The  equations  of  mo¬ 
tion  of  the  object  are  then  as  follows  from  application  of  Newton's 
second  law  to  the  problem  (e.g.,  refs.  1  and  2): 


^  A 

dt2  r^  "  ^ 


A 

dt^  r^  y 
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(1) 
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In  the  notation  adopted  hereinafter,  equation  (l)  will  he  indicated  as 


X  y,z 


(1) 


where  x  y,z  Indicates  that  separate  equations  in  x,  y,  and  z  are 
included  to  complete  the  set.  All  symbols  are  defined  in  appendix  A. 

The  acceleration  term  p.x/r^  is  due  to  central  gravitational  attraction, 
where 

M.  =  k^(MQ  +  m)  (2) 

(k^  is  the  gravitational  constant,  is  the  mass  of  the  body  at  the 

origin,  and  m  is  the  object's  mass).  The  components  of  the  disturbing 
acceleration  that  disturb  the  two-body  orbit  are  A^,  x  ->  y,z  in  equa¬ 
tions  (l)  and  may  be  written  as 


n 


i=l 


where  is  the  mass  of  the  i'^h  disturbing  body,  n  is  the  number  of 

gravitating  bodies  excluding  the  central  body,  is  defined  by 

At  :=  (x  -  x^)^  +  (y  -  y^)2  +  (z  -  (4) 

and  f^  is  the  component  of  the  disturbing  acceleration  along  OX  due 

to  all  other  forces.  For  example,  these  may  include  propulsion  thrust, 
aerodynamic  forces,  and  forces  due  to  the  oblateness  of  M^.  Wo  restric¬ 
tion  need  be  placed  on  the  form  of  disturbing  acceleration  except  that 
it  be  sufficiently  well  defined  to  permit  integration. 

Equations  (l)  may  be  integrated  directly  by  numerical  methods.  How¬ 
ever,  in  many  cases  larger  intervals  may  be  used  in  numerical  integra¬ 
tion,  or  approximate  closed-form  solutions  can  be  obtained  if  the  equa¬ 
tions  are  expressed  in  terms  of  perturbations  of  orbital  elements. 


The  perturbation  theory  uses  as  a  reference  an  orbit  having  no  per¬ 
turbations.  If  the  disturbing  acceleration  is  assumed  to  be  zero,  the 
differential  equations  become 


d  X 

dt^  r^ 


0, 


X  ->  y,z 


(5) 
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The  solution  of  eq^uations  (s)  is  readily  obtained  and  is  found  to  be  a 
conic  section.  The  motion  of  the  object  can  be  represented  by  six  or¬ 
bital  elements  obtained  from  the  constants  of  integration  of  equations 
(5) .  Table  I  gives  a  collection  of  two-body  equations  relating  selected 
orbital  elements  and  parameters.  Although  it  is  not  possible  to  express 
the  Cartesian  positions  and  velocities  explicitly  in  terms  of  the  orbital 
elements^  the  solutions  of  equations  (s)  may  be  indicated  as 


where 


X  =  x(c3^,  Cg,  C3,  C4,  Cg,  Cg,  t). 

x  y,z 

(6a) 

X  —  x(c4,  Cg,  Cg,  C4,  Cg,  Cg,  t). 

X  ^  y,z 

(6b) 

Cg,  •  •  ‘ are  orbital  elements. 


If  A^  0,  X  y, z  in  equations  (l),  the  path  is  not  a  simple 

conic.  However j  at  any  instant  it  may  be  regarded  as  a  conic  with  vari¬ 
able  orbital  elements.  In  fact,  equations  (6)  are  the  solutions  of  equa¬ 
tions  (1)  if  the  orbital  elements  c^^,  c^,  •  •  Cg  are  regarded  as 

variables . 


This  introduces  the  concept  of  the  oscvaating  orbit.  Let  an  object 
be  moving  in  a  perturbed  path  about  a  central  body.  An  instantaneous 
two-body  orbit  always  exists  tangent  to  the  actual  path  at  the  point  and 
having  a  velocity  in  the  orbit  equal  to  that  of  the  actual  body.  Such 
tangent  orbits  are  called  the  osculating  orbits.  The  relations  implicit 
in  this  definition  are  used  to  derive  the  equations  for  the  disturbed 
orbit  in  terms  of  orbital  elements.  The  three  second-order  differential 
equations  of  equations  (l)  are  transformed  to  six  first-order  simultaneous 
differential  equations  involving  orbital  elements.  Lagrangian  brackets 
are  utilized  to  solve  the  set  of  simultaneous  differential  equations  for 
derivatives  of  the  various  orbital  elements  explicitly.  The  explicit 
derivatives  are  the  objective  of  this  report. 


As  the  first  step  in  obtaining  these  derivatives,  equations  (6)  are 
differentiated  regarding  0^,0^,  •  •  •,  Cg  as  variables.  The  following 

equations  that  apply  to  the  actual  path  are  obtained: 


dx  ,  \  5x 

dt  =  5t  ^ 


/  v^^k 
k=l 


dci 


Scv  dt  ' 


(7a) 


d^x  _  Sx  _j,\  '  dx  ^ 
dt2  ^ 

k?=l 


X  -»  y,z 


(7b) 
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Ec[uations  (l)  m^y  be  introduced  to  eliminate  (d^x/dt^),  x  -»•  y,z  from 
eq_uations  (7b).  Then  eq_uations  (?)  may  be  written  as 
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Is  _  ^  ^  =  0 

§t^  dt  J  ’ 


X  -»  y,z 


6 


k=l 


dt 


X  -►  y,z 


(8a) 


(8b) 


In  the  osculating  orbit  dcj^/dt  =  0  and  -  0,  so  that  equations 
(S)  become 

S¥'dt  =  °'  X  -  y,z  (9a) 

||  +  ^=0,  x-y,z  (9b) 

Introducing  the  requirement  that  velocities  in  the  actual  path  and 
in  the  osculating  orbit  are  equal,  equations  (9a)  may  be  substituted  into 
equations  (8a) j  and  similarly,  because  the  acceleration  in  the  osculating 
orbit  differs  from  that  of  the  true  orbit  only  by  the  disturbing  acceler¬ 
ation,  equations  (9b)  may  be  substituted  into  equation  (8b)  to  give 


(lOa) 


X  -►  y,z 


(10b) 


Equations  (lO)  are  the  resulting  six  first-order  differential  equations. 
They  are  not  ad.apted  for  computation  because  equations  (6)  are  not  explic¬ 
itly  available  and  because  the  derivatives  of  the  orbital  elements  appear 
simultaneously  rather  than  explicitly.  This  difficulty  is  conveniently 
removed  by  further  manipulation.  The  following  equations  are  written  in 
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a  form  convenient  for  formulation  in  terms  of  Lagrangian  brackets: 


(11a) 


The  validity  of  equations  (lla)  is  obvious  because  each  term  contains  a 
zero  factor  from  equations  (lO).  Results  of  the  operations  shown  may  be 
written  as 


^  ~  ^c^^  ^  ~  1>2^’  •  ',6 


(lit) 


k=l 


where 


[c.,c.]  =  +  + 

J  dc.  dcn,  dc.  dc-^  dc^i  dc^  dc-  dc^  dc-  dcv  ocT  o 


dcv 


and 


I>Cj  =  If-  +  Ay  ll:  +  A,  If-,  J=l,2,- 

j  j  j 

The  brackets  [c^Cj^]  are  the  Lagrangian  brackets . 


:i2) 


(13) 
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General  Formula  for  a  Lagrangian  Bracket 


Previously^  a  general  set  of  orlital  elements  has  been  used.  Before 
proceeding  to  evaluate  the  Lagrangian  Brackets  of  equations  (llh),  it  is 
convenient  to  choose  a  specific  set  of  independent  orbital  elements  so  as 
not  to  complicate  the  analysis.  The  set  chosen  is  not  significant;  for, 
as  will  he  shown  in  the  RESULTS  section,  it  is  relatively  simple  to  sub¬ 
stitute  any  elements  that  may  be  desired.  If  there  are  chosen  as  the  set 
of  orbital  elements  the  semimajor  axis  a,  eccentricity  e,  time  of  peri- 
center  passage  t^,  argument  of  pericenter  05,  orbit  inclination  I,  and 

longitude  of  the  ascending  node  Q,  the  expression  for  the  Lagrangian 
bracket  is 


[s,q] 


\  "'P'*  2a  )  ,  (o5, -\/M-a(l  -  e^)  )  ,  (fi,-»/pa(l  -  e^)  cos  l) 
(s,a)  *  (s;^) - ^  (s^,) - 


(14) 


where  s  and  q  are  any  of  the  orbital  elements.  The  right-hand  side 
is  expressed  in  Jacobian  notation.  Note  that  [s,s]  =  0  and 
[s,q]  =  -[q,s]  from  equation  (12) .  The  derivation  of  equation  (14)  is 
given  in  appendix  B.  It  results  from  geometric  relations  existing  among 
the  instantaneous  values  of  the  orbital  elements  of  any  orbit.  It  will 
be  used  to  evaluate  the  Lagrangian  brackets  of  equations  (11b) . 


Evaluation  of  Lagrangian  Brackets 

A  Lagrangian  bracket  appears  in  equations  (lib)  for  each  of  the  36 
combinations  of  the  six  chosen  orbital  elements.  Evaluation  from  equa¬ 
tion  (14)  shows  six  of  the  twelve  nonzero  brackets  to  be 


[a,tp]  =  - 


a(-p/2a) 
^t^^ - 


(15a) 


[a,J2] 


[a, 05]  =  -  ^VH.a(l,.-_e21  , 

doa  da 


dSi  d-/|ia(l  -  e2)cos  I 
Sii  Sa 


(15b) 


(15c) 


doo  d-/ua(l 
.  de 


[6,05]  -= 


(I5d) 
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[e,S2] 


^  ^-ZliaCl  -  e^)~cos  I 

Sir  §e 


e-y/^a 

-y/l  - 


cos 


I 


[I;^] 


SQ  S-y/p.a(l  -  e^~cos  I 

Sir  5i 


Vl^a(l  -  e2) 


sin  I 


(I5e) 


(I5f) 


By  observing  the  property  of  the  brackets  that  [s.q]  =  -[q,s],  the  values 
oi  the  remaining  six  nonzero  brackets  are  apparent  from  equations  (l5). 
Omitting  all  zero  brackets,  equations  (lib)  become 


[a,tp]tp  +  [a^cojclj  +  [a,S2]d  =  (l6a) 

[e,ai]<i)  +  (let) 

=  Dj  (l6c) 

-[a,tp]a  =  (l6d) 

-[a,C0]a  -  [e,CD]4  =  (l6e) 

-[a,fl]a  -  [e,S2]4  -  [l,S2]i  =  (l6f) 


Derivatives  of  Elements  in  Terms  of  Disturbing 
Functions  of  Elements 

Introducing  the  values  obtained  for  the  brackets  in  equations  (15) 
into  equations,  (l6)  and  solving  the  system  for  the  derivatives  give  the 
following  form  for  the  Lagrangian  equations : 


a  = 


a(l  -  e^) 
pe 


(17a) 


e 


(17b) 
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I  =  '  /  .  .=  cot  I  -  —  _ 

Vl^ad  -  e2)  \  “  sin  I 


-f,  =  ~  ^  ^  D 

i-ie  e  ^ 


VI  -  e2  ^ 

-  ^(= 

|j.a  s 


cot  I 


■V|-ia(l  -  e2) 


D-, 


Q.  = 


sin  lVpa(l  -  e2) 


(17c) 


(I7d) 


(I7e) 


(I7f) 


Disturbing  Functions  of  Elements  in  Terms  of 


Components  of  Disturbing  Acceleration 

Prior  to  the  integration  of  equations  (17)  it  is  desirable  to  express 
the  disturbing  functions  Dg^,  Dg^  D^  ,  D^,  Dj^,  and  Dj  in  more  convenient 

P 

form. 


If  s  is  any  of  the  elements  a,  e,  tp^  co,  or  I,  any  of  equa¬ 
tions  (13)  may  be  written  as 


D. 


(13a) 


It  is  necessary  only  to  evaluate  dx/Ss^  x  ->  y, z  for  each  element  in 
order  to  obtain  useful  forms  of  equations  (l7)  that  may  be  integrated 
(either  formally  or  numerically).  However,  the  equations  will  reduce  to 
a  more  convenient  form  if  the  Cartesian  disturbing  acceleration  components 
A^;  Ay,  are  resolved  into  a  new  orthogonal  set  (fig.  2)  as  follows: 

(l)  a  component  normal  to  the  orbital  plane  W,  positive  when  A^  is 

positive;  (2)  a  component  normal  to  the  radius  and  in  the  orbital  plane 
C,  positive  when  making  an  angle  of  less  than  n/z  with  the  direction 
of  motion;  (3)  a  component  along  the  radius  R,  positive  when  pointing 
outward  from  the  origin. 


The  analysis  in  appendix  C  determines  the  disturbing  functions  Dg 
Dg,  D^  ,  D^,  0^2,  and  Dj  in  terms  VJ ,  C,  and  R  as  shown  in  equations 

(18)  .  ^ 


■p  IP 
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T  r  cos  u  ,, 
1  =  7=z= —  w 


(I9f) 


where  p  _  is  'che  semilatus  rectum.  These  results,  together  with  other 
lorms  derived  from  them,  are  given  in  table  II. 

Alternate  Components  of  Disturbing  Acceleration 

Components  of  the  disturbing  acceleration  in  the  orbital  plane  mav 
be  alternately  taken  tangent  to  the  path  T  and  normal  to  the  path  W; 

IS  positive  in  the  direction  of  motion,  and  N  is  positive  when  di- 
rected  toward  the  interior  of  the  orbit.  Substitution  in  terms  of  T 

^  ^  introducing  expressions  from  table  III  gives 

tne  following  changes  in  the  derivatives  of  equations  (l9) : 


p  V 

a  =  2a^  -  T 


VL  e 


e  =  U2(e  +  cos  v)T  -  jl-Cl.  "  e^)sln  v  . 


(20a) 


(20b) 


sin  V 


T  +  (2e  +  cos  v  +  e^  cos  v)W  -  (-^= 

^  VVu.. 


sin  u  cot  I 


(20c) 


(e^  +  e  cos  V  +  l)sin 


-  3(t  -  tp)(l  +  e2  +  2e  cos  v)  ^  +  (f  cos  v^nJ 


(20d) 


where  V  is  the  velocity. 


Elimination  of  (t  -  tp) 

The  quantity  (t  -  tp)  may  be  eliminated  from  equation  (20d)  by  in¬ 
troducing  table  I  equations  (l-96)  and  (l-97)  for  the  cases  e  <  1  and 
e  >  1,  respectively.  The  results  for  the  two  cases  are 


€  L  ^  e2 


sin  V 


3(1  +  2e  cos  V 

Vl  -  e2 


+  (^  cos  v)tlj, 


0  <  e  <  1 
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where  E  is  the  eccentric  anomaly,  anr^ 

-  tT  f +  3e  +  +  cos  v)! 

P  ^>Pyi-e2^|_e  1  +  e  cos  v  J 


Sin  V 


5(1  +  2e  cos  V  +  e^) 


■7 


e^  -  1 


e  >  1 


(22) 


where  F  -  -iE  in  the  hyperbolic  orbit  corresponds  to  E  in  the  elliptic 
orbit.  Thus,  F  is  imaginary  when  an  elliptic  orbit  exists,  and  E  is 
imaginary  in  the  case  of  the  hyperbolic  orbit.  The  result  for  e  >  1  is 
identical  to  the  result  for  e  <  1,  but  with  E  replaced  by  iF  and 
-y/l  _  replaced  by  i-/e2  -  l. 


The  value  for  tp  when  e  =  1  is  not  directly  evident  as  the 
quantity  (l  -  e  )  ^  0  when  e  1.  Equation  (l9c)  for  t^  may  be 
written  as  ^ 


= 


(1  -  e^)  [  Ll  +  e  cos  V 


I  sin  V  2  + 
[  e  FT" 


e  cos  V 


e  cos  V 


cos  V 


3(1  + 


3e  sin  V  (t  -  tp)J  R 


e  cos  V, 


(t  -  tp)|  C 


(23) 


Eliminating  (t  -  tp)  from  equation  (23)  with  table  I  equation  (1-95) 


yields 


P 


■  1.(1  -  e2: 


cos  V 


3e^sln^v 


1  +  e  cos  V  e  7t  2\  !  T 

(1  -  e^) (1  +  e  cos  v) 


3e  sin  v 
1  -  e2 


f 

/-V 


dv 


1+e  cos  V 


R  + 


sin  V  / 2  +  e  cos  v\ 
Vl  +  e  cos  v) 


+  5e  sin  v  3(l  +  e  cos  v) 
1  -  e2  1  -  e2 


/ 


dv 


1+e  cos  V 


(24) 
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Using  the  relation 


•  V 


cos  V  dv 


(l  +  e  cos  v)2  2(1  -  e^) 


sin  V 


(1  +  e  cos  v)  ^ 


+  (1  t  2e^)sin 


3e 


(1  -  e^)(l  +  e  cos  v)  1  -  e^  /  1  +  e  cos  v 


r 


dv 


to  eliminate 


dv 


e  cos  V 


Jc 


from  equation  (24)  yields 


P  |i 


cos  V  dv 


cos  V 


(1  +  e  cos  v)^  e(l  +  e  cos  v)^ 


+  I  —  (1  +  e  cos  v) 


/ 


cos  V  dv 


(1  +  e  cos  v)^ 


Equation  (26)  is  now  defined  at  e  =  1,  since 


f 

J  n 


cos  V  dv  _  sin  v  1+5  cos  v  +  cos^v 
(1  +  cos  v)2  5 


(1  +  cos  v)^ 


Hence^  the  equation  for  t  on  a  parabola  is 


t„  = 


p  5,(r/c03  v)2  { [2(1  +  =  “O'*  -  +  =1"  -]  C 

+  js  -  cos  V  -  4  cos^v  -  2  cos^-\^ 


(25) 


(26) 


(27) 
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Derivatives  of  Alternate  Orbital  Parameters 

to 

The  element  "p"  may  be  used  in  place  of  the  element  "a"  as  a  variable 
of  integration.  Similarly^  the  pericenter  radius  rp  may  replace  a  or 

p.  As  p^  a.,  and  rp  are  similar  elements,  the  q^uestion  of  which  to 
select  is  determined  by  numerical  and  convenience  considerations. 

Semilatus  rectum.  -  Taking  the  derivative  of  p  =  a(l  -  e^)  yields 

p  =  a(l  -  e^)  -  2eae.  Introducing  expressions  for  a  and  e  from  eaua-  W 

tions  (19a)  and  (l9b)  yields 

p  =  2r  C 

T  p 

or,  in  terms  of  T  and  K, 

p  =  ^  e(sln  v)l^ 

.  -  The  equation  for  radius  of  pericenter  is 
=  p/(l  +  e) ,  and  its  derivative  is 

-  "■p®)  (30) 

Introducing  expressions  for  p  and  e  from  equations  (29)  and  (2ab)  into 
equation  (30)  and  simplifying  yield 

"  7[^^p  ^  1  +°e  ^  +  (^  sin  v)n]  (31) 

Mean  anomaly.  -  Also,  it  may  be  desirable  to  use  the  mean  anomaly  M 
in  place  of  the  element  tp.  Here,  the  distinction  is  significant  in 

that  M  is  not  an  element  and  varies  even  on  a  two-body  orbit.  From  the 

familiar  expressions  M  =  n(t  -  tp)  and  n  =  p/ a^  in  table  I,  the 

derivatives  become 

Si  =  n(t  -  tp)  +  n  -  ntp  (32) 

and 


(28) 


(29) 


© 


202T-: 
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where  equation  (20a)  was  used  for  a  in  equation  (33) .  Now^ 


(34) 


when  tp  is  taken  from  equation  (20d). 

Introducing  equations  (33)  and  (34)  into  equation  (32)  yields 

Vl  -  lo/  •  ,  l\m  r  cos  V  „ 

M  =  n  -  ^ -  |2  (sin  v)^—  +  -jl  +  -  nJ  (35) 

It  should  he  noted  that  M  and  M  are  both  zero  when  e  =  1  and  that 
M  reduces  to  its  two-body  value  when  the  in-plane  components  of  the 
disturbing  acceleration  are  zero. 

True  anomaly .  -  The  orbital  parameter  y,  true  anomaly,  may  be  used 
as  an  alternate  to  either  the  element  tp  or  the  parameter  M.  From 

table  I  equations  (1-47),  (1-18),  and  (1-24),  Kepler's  equation  may  be 
written  as 


,,  ^  -1  Vl  -  e2  sin  V  -i/i  -  e^  sin  v 

M  =  tan  — - - e  V-r - 

e  +  cos  V  1  +  e  cos  v 

Taking  derivatives  in  the  disturbed  orbit  yields 


(36) 


M  =  — - r  rv(l  -  e^)  -  e  (sin  v)(2  +  e  cos  v)"]  (37) 

(1  +  e  cos  v)^  L  J 

Solving  equation  (37)  for  v  and  introducing  equations  (35)  and  (20b) 
to  eliminate  M  and  e,  respectively,  yield  the  result 


2e  +  e^  cos  v  +  cos  v 
Ve(l  +  e  cos  v) 


(38) 


As  in  the  case  of  M  (eq.  (35)),  it  will  be  noted  that  equation  (38)  re¬ 
duces  to  the  two-body  derivative  expression  when  the  in-plane  components 
of  the  perturbative  acceleration  are  zero. 
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Table  IV  Is  a  q^ualitative  summary  showing  how  the  components  of  the 
disturbing  acceleration  affect  the  derivatives  of  the  various  orbital 
elements  and  parameters. 


Orbital  Element  Relations  for  a  Circular  Orbit 


In  a  circular  orbit  the  location  of  the  pericenter  is  undefined. 
Conseq^uently ^  the  elements  cd  and  tp  and  the  true  and  mean  anomalies 

V  and  which  are  related  to  the  location  of  the  pericenter,  are  un¬ 
defined.  Thus,  relations  involving  a>,  tp,  v,  and  M  take  on  an  in¬ 
determinate  form.  However,  any  perturbation  having  a  component  in  the 
plane  of  the  orbit  will  immediately  establish  the  limiting  values  of  co. 


tp,  V,  and  M  in  the  circular  orbit 


The  circular  orbit  expressions  are 

derived  for  the  T,  E  system  of  resolution  of  the  in-plane  disturbing 
acceleration  components .  The  derivatives  are  also  given  in  terms  of  the 
C,  R  system  in  table  II. 


True  anomaly.  -  Taking  the  table  I  equation  (1-117)  for  cos  v, 

e  =  0 


(p/r)  -  1 

cos  V  =  — - - 


and  applying  L ‘Hospital's  rule  give 

(rp  -  pr)/r^ 

cos  V  =  - , 

e 


e  =  0 


(39) 


Substituting  for  p  from  equation  (29),  for  e  from  equation  (20b), 
and  for  f  from  table  I  equation  (1-80)  [eq.  (1-80)  is  also  valid  for 
disturbed  orbits]  into  equation  (39)  yields  in  the  limit  as  e  approaches 
zero 


cos  V 


_ 2T _ 

2T  cos  V  -  N  sin  v’ 


e  =  0 


(40) 


Equation  (40)  is  valid  if  e  =  0  and  e  ^  0.  The  latter  condition  re¬ 
quires  that  an  in-plane  ccmponent  of  the  perturbing  acceleration  exists. 
Solving  for  sin  v  and  cos  v  from  equation  (40)  yields 


N 

sin  V  =  +  — ’  ■  — . 

Y4T2  +  n2 


e  =  0 


(41a) 


cos  V 


2T 


_y4T2  +  k2 


e  =  0 


(41b) 

*  X; 
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where  the  upper  sign  gives  the  limiting  value  of  the  true  anomaly  for 
leaving  a  circular  orhit  and  the  lower  sign  gives  the  limiting  value 
when  entering  a  circular  orhit.  It  can  he  seen  from  equations  (4l)  that, 
if  a  perturbing  force  is  arranged  to  force  an  initially  elliptical  orhit 
through  circular,  the  value  of  true  anomaly  will  make  a  step  change  of 
Tt  radians . 

Argument  of  pericenter.  -  The  argument  of  perlcenter  is  determined 
from  table  I  equation  ( 1-176) : 

CD  =  u  -  V 

Time  of  pericenter  passage.  -  In  a  circular  orhit  the  time  from 
"pericenter,"  t  -  tp,  equals  the  arc  length  to  "pericenter"  divided  by 

the  velocity;  that  is. 


t  -  tp  =  ^  V, 


e  =  0 


Derivatives  of  Orbital  Parameters  in  a  Circular  Orbit 

Semllatus  rectum  and  semima,ior  axis.  -  Equations  (29)  and  (20a)  for 
p  and  a  reduce  directly  to  the  same  expression  when  e  =  0: 

P=4=2:^T,  e=0  (43) 

Eccentricity  and  radius  of  pericenter.  -  Using  the  limiting  values 
for  sin  v  and  cos  v  from  equations  (41)  in  equations  (20b)  and  (3l) 
for  e  and  r  yields  the  following: 

hr 

e  =  ±  i  V4:T2  +  ii2^  e  =  0  (44) 


r  =  y  (2T  +  V4T2  +  N2)^  e  =  0 

Jr  V 


where  again  the  upper  signs  give  the  limiting  values  of  e  and  fp  when 
leaving  a  circular  orbit  and  the  lower  signs  are  for  entering  a  circular 


orbit . 


Argcmient  of  pericenter.  -  Equation  (20c)  may  be  rewritten  to  give 
CD  =  ip-(2T  sin  V  +  -  K  cos  v)  +  -  (2  +  e  cos  v)N  -  (  sin  u  cot  I|W 

P  /  p  VVt^p  J  . 


(46) 
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and  tp  for  a  circular  orbit  reduces  to 


;  _  1  ,  r/TK  -  NT 


r  /ST^  +  N^ 


2  \  /im2 


P  2  V^^2  +  +  u 

Mean  anomaly.  -  Equation  (35)  may  be  written  as 


N  +  3vT  L  e  =  0 


(51) 


M  =  n 


Vl  -  e2 

V 


2(sln  v)T  +  —  (l  -  e^)N  cos 


V 


+  2 


(*F  '■)' 


(52) 


The  limi of  the  indeterminate  term  in  equation  (52)  is 


Limit  2T  sin  V  +  I  (1  -  e^)N  cos  v  _  ^  ^  ^^2 

e  -*  0  e  ----  -  N 


2TW 


2r  \4^2  4  jj2y  ^41.2  +  n2 

TN  -  NT 


+  V 


4T2  +  n2, 


and  M  for  a  circular  orbit  becomes 


I_  +  br/6T2  +  n2^  ^  2TN  TN  -  NT 

2r  V\4t2  +  ^2j  "  +  ig2  '  41^  +  n2 


(53) 


where  the  upper  sign  is  for  leaving  a  circular  orbit  and  the  lower  sign 
is  for  entering  a  circular  orbit. 

True  anomaly.  -  Equation  (38)  may  be  written  as 


1  fl  /, 

-  V 


2T  sin  V  +  —  N  cos  v 
P  / 


—  (2  +  e  cos  v)N  + 

Vp  "  r 


Vpp 

2 


which  yields  in  the  limit  as  e  approaches  zero 


V  = 


N/2T^  +  N^ 


2r  ^\4'ii2  4  ]\f2 


TN  -  NT 
4T^  4 


e  =  0 


(54) 
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A.^plication  of  derivatives  for  circular  orbits.  -  If  the  disturbing 
force  is  applied  in  a  constant  direction  relative  to  the  velocity  direc¬ 
tion,  it  may  he  noted  that 

TN  -  NT  =  0 

even  if  the  magnitude  of  the  disturbing  acceleration  is  varying.  Thus, 
if  the  direction  of  the  disturbing  acceleration  is  fixed  relative  to  the  t-J 

velocity,  all  terms  involving  N  and  T  will  vanish  in  equations  (49)  ^ 

to  (54) .  Similarly  if  the  direction  of  the  disturbing  acceleration  is  o 

fixed  relative  to  the  radial  direction,  all  terms  Involving  C  and  R 
v;ill  vanish  from  the  equations  in  table  II. 

The  equations  for  circular  orbits  given  in  the  preceding  sections 
and  in  table  II  are  completely  valid  only  for  circular  orbits.  However, 
certain  of  the  equations  will  be  found  to  be  sufficiently  accurate  for 
near-circular  orbits  and  can  be  used  as  the  basis  of  approximate  equa¬ 
tions  . 


CONCLUDING  REMARKS 

All  the  results  derived  herein  for  the  perturbation  derivatives  of 
the  various  orbital  elements  and  parameters  are  listed  in  table  II. 

Table  III  contains  expressions  for  the  orthogonal  system  of  components 
C,  R,  and  VI  in  terms  of  the  Cartesian  components  of  the  disturbing 
acceleration  A^,  A^.,  A^  in  the  OX,  OY,  OZ  system.  It  also  contains 

equations  interrelating  the  C,  R  and  T,  N  .'systems  of  the  in-plane 
components  of  the  disturbing  acceleration.  Table  I  is  a  collection  of 
various  forms  of  the  two-body  equations  that  also  apply  to  the  osculating 
orbit . 

As  illustrated  for  p,  fp,  M,  and  v,  perturbation  derivatives  of 

other  alternate  elements  and  parameters  may  be  derived  from  the  expres¬ 
sions  in  tables  I,  II,  and  III.  Other  integration  variables  that  nay 
be  useful  for  the  solution  of  problems  in  orbital  mechanics  are  suggested 
in  references  2,  3,  and  4. 

Selection  of  the  best  set  of  orbital  elements  or  parameters  for  a 
particular  type  of  special  perturbation  problem  in  orbital  mechanics  de¬ 
pends  on  the  nature  of  the  problem.  However,  it  is  expected  that  examina¬ 
tion  of  the  derivatives  given  herein  will  help  to  Indicate  which  param¬ 
eters  should  be  used  for  a  specific  problem. 
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APPENDIX  A 

SYMBOLS 


The  following  symbols  are  used  in  this  report; 
component  of  disturbing  acceleration 

a  semimajor  axis  of  conic  section,  negative  in  hyperbolic  case 

b  semiminor  axis  of  conic  section 

0  perturbative  acceleration  in  circumferential  direction 

c  constant  of  integration 

disturbing  acceleration  function  for  element  s, 

.  Sx  ,  .  3y  ,  .  Sz 


E 

Eg 

e 

F 


G 

E 


eccentric  anomaly 
energy  per  unit  mass 
eccentricity 

used  in  hyperbolic  orbits  to  correspond  to  eccentric  anomaly 
in  elliptic  orbits,  F  =  -iE 

component  term  of  disturbing  acceleration  includes  forces 
due  to  all  except  gravitating  bodies  not  located  at  prob¬ 
lem  origin 


a2/2(E  -  e  sin  E) 

-  lil/2al/2(3E  +  e  sin  E) 


h 

I 

i 


angular  momentum  per  unit  mass  equals  twice  the  rate  of 
description  of  area  in  orbital  plane 

orbital  plane  inclination 

unit  imaginary  number,  n/ -1 


X 


5X 

6s 


+  Y 


6y 

55 
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M 


Mi 


m 


N 


n 

P 

P 


Q,S 


q,s 

R 


Rg 


r 


r . 

1 


T 

t 

U,Ui 

u 


gravitational  constant 
mean  anomaly 
gravitating  body  mass 

gravitating  body  mass  at  problem  origin 
object  mass 

perturbative  acceleration  in  orbital  plane  in  direction 

normal  to  velocity,  positive  when  directed  toward  interior 
of  orbit 

mean  angular  orbital  motion  of  object,  2n/P 

orbital  period 

semilatus  rectum,  a(l  -  e^) 

denote  functions  defined  for  convenience 

any  pair  of  orbital  elements 

perturbative  acceleration  in  radial  direction,  positive 
outward 

range  on  surface  of  sphere  intersected  by  an  elliptical 
orbit 

radius  from  origin  to  object 
radius  from  origin  to  disturbing  body 

pericenter  radius,  the  minimum  distance  from  central  body 
to  orbit 

radius  of  sphere 

perturbative  acceleration  in  direction  of  velocity 
time 

time  of  pericenter 

complex  variables,  x  +  iy,  x^  +  ly^ 

argument  of  latitude,  angle  measured  from  ascending  node  to 
object  radius  in  the  direction  of  motion,  u  =  v  +  oo 
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P]_ 


'^2^P2jTf2 


object  velocity  or  complex  variable^  x  -  iy 
complex  variable^  x^^  - 

true  anomaly^  angle  measured  from  perlcenter  to  object  radius 
in  the  direction  of  motion 

perturbative  acceleration  normal  to  orbital  plane ^  positive 
in  the  direction  of  axis  OZ 

denote  object  coordinates  in  the  OYg  axis  system 

vhere  OXj  coincides  with  pericenter  radius  and  OY^  lies 
in  orbital  plane  (see  fig.  2) 

denote  object  coordinates  in  the  OX^  OY^  OZ  Cartesian 
system  (see  fig.  l) 

direction  cosines  of  the  OX3  axis  referred  to  the  OX^  OY^ 
and  OZ  axes,  respectively 

direction  cosines  of  the  OY3  axes  referred  to  the  OX,  OY, 
and  OZ  axes,  respectively 

distance  from  object  to  any  pertxirbing  gravitating  body 
natirral  logarithm  base 

gravitational  constant  equals  acceleration  of  object  at  unit 

distance  from  M  due  to  M  ,  u.  =  k2(M  +  m) 
o  o  o 

used  in  hyperbolic  orbits  to  correspond  to  n  in  elliptic 
orbits,  V  =  -in 

path  angle,  angle  between  circumferential  and  velocity  di¬ 
rections,  positive  clocta>rise 

arg;rment  of  pericenter,  angle  measured  from  ascending  node 
to  pericenter  radius  in  direction  of  motion 
longitude  of  ascending  node 


Subscripts: 


apocenter  conditions  in  table  Ij  elsewhere,  disturbing  func¬ 
tion  of  semimajor  axis 

disturbing  function  of  e 
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I  disturbing  function  of  I 

i=l,2,'-',n  disturbing  gravitating  body  number 
indicates  equation  of  a  set 

indicates  orbital  element 

p  pericenter  conditions 

s  conditions  of  spherical  surface 

tp  distirrbing  function  of  t^ 

x-y,z  indicates  extension  of  equations  in  x  to  a  system  in¬ 

cluding  y  and  z 

x,yjZ  components  taken  about  or  along  an  axis 

CO  disturbing  function  of  co 

1,2,3  coordinate  system  and  coordinates  in  them 

disturbing  function  of  £2 

Superscript : 

•  indicates  derivative  with  respect  to  time 
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WHITTAKER'S  DERIVATION  OF  THE  GENERAL  FORMULA 
FOR  A  LAGRANGIAN  BRACKET 

The  following  development  is  applicable  to  the  equations  of  any 
osculating  orbit.  It  is  similar  to  that  given  in  article  5-16  of  Smart 
(ref.  1) . 

Let  s  and  q  be  any  two  of  the  six  orbital  elements  and  x^  y,  z 
be  the  object's  coordinates;  then  equation  (12)  may  be  written  as  a  sum 
of  Jacobi ansj 


[s,q] 


(s,q)  (s,q) 


(z,z) 

( s,q) 


(Bl) 


where j  by  definition^ 


(x,x)  _  dx  Sx  5x  3x 
( s,q)  ^  ^  Sq 

and  so  forth,  for  the  other  Jacobians.  The  objective  is  to  transform 
equation  (Bl)  into  an  expression  of  the  same  form  but  involving  only  ex¬ 
plicit  functions  of  the  orbital  elements.  Define 


and 


Then, 


(B2) 


[s,q] 


5q  ds 


(B3) 


Let  the  OX,  OY,  OZ  axis  system  in  figure  1  be  fixed  to  some  reference 
astronomical  line  and  plane  such  as  the  mean  equinox  and  equator  of 
1950.0.  The  orbit  plane  is  projected  on  the  celestial  sphere.  The 
equatorial  longitude  of  the  ascending  equatorial  node  referred  to  the 
vernal  equinox  is  then  denoted  by  S,  and  the  inclination  of  the  orbital 
plane  to  the  equator  is  I.  Rotate  OX  and  CY  about  OZ  through  the 
angle  n  to  obtain  the  axis  system  OX^,  OY^,  OZ^ .  Coordinates  of  the 
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object  in  the  OX^  OY^  OZ  system  (x^y^z)  are  given  in  terms  of  those  in 
the  OZ^  system  (x^,y^,Zj^)  by 


X  =  x^  cos  -  y^  sin  ^2 
y  =  X j_  sin  S2  +  y^  cos  SI 


It  is  convenient  to  introduce  complex  variables  by  the  following  defini¬ 
tions  ; 


U=x+iy^  V=x-iy 


Ui  =  +  iyj_,  =  x^  -  iy^  J 

where  ^  i  =  Then^  by  using  equations  (B4)  with  the  definitions  in 

equations  (B5),  it  may  be  shown  that 


U  =  U2^(cos  S2  +  i  sin  S2) 
V  =  Vj_(cos  S2  -  1  sin  S2) 


(B6a) 

(B6b) 


U  =  Uqe^ 


(B7a) 


V  =  V  p' 


(B7b) 


The  time  derivative  of  equation  (B7a)  and  the  partial  derivative  of  equa¬ 
tion  (B7b)  with  respect  to  s  are^  respectively. 


U  =  Uqe^ 


Ss  ^  ( ^s  "  ^  1  ds 


Multiply  equations  (B8)  together  to  obtain 


^  d^  -  ^1  5^  -  ^^1^1  51 
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Taking  the  real  parts  of  equation  (B9)  yields 


(xpyi  -  Xiy3_)  ^ 


(BIO) 


But  it  is  seen  from  equations  (l-33)  and  (l-44)  of  the  two-hody  orbital 
relations  in  table  I  that  x^yj^  -  x^y^^  =  h  cos  I,  which  is  the  projection 

on  the  -  OYj_  plane  of  twice  the  rate  of  description  of  area  in  the 

orbit  plane.  Hence,  from  equations  (B2)  and  (BIO) 


Sx]  Sy-i  Sz-i  ;5iq 

S  =  ^1  ^  +  yi  ^  +  h  cos  I  g  (Bll) 

as  z  =  from  equation  (B4)  .  Now,  rotate  the  OY^,  -  OZ^  plane  about 
OXp  through  the  angle  I  to  obtain  the  axis  system  0X3,  0X3,  OZg  hav¬ 
ing  object  coordinates  (x2,y2^Z2),  so  that  the  OX2  -  OY2  plane  coincides 
with  the  orbit  plane.  Then,  by  analogy  with  equation  (BIO), 


.  Syi 
^1  ST- 


dz^ 
"1  ds 


.  ay2  3z2  .  xj 

^2  57"  ^2  ^  (^2^2  -  ^2^2^  Si 


(B12) 


But,  since  X2_  =  X2  and  Z2  =  0,  equations  (Bll)  and  (B12)  nay  be  used 
to  yield 

.  5x2  Sy2  ^ 

s  =  ^2  ST"  ^2  ST"  ^  ^  Si 

Now  rotate  the  OX2  -  OY2  plane  in  the  orbital  plane  through  the  angle 
03  to  obtain  the  axis  system  OX7,  OY3  having  object  coordinates  (X,Y) 
so  that  0X3  lies  on  the  pericenter.  Again,  by  analogy  with  equation 
(BIO)  it  follows  that 


Sxc 


^y2  sx 


Sy 


Sco 


^2Sr"y2Sr=X^-tYgtH-(XY-YX)^ 

But  XY  -  YX  =  h;  thus,  use  of  equations  (B13)  and  (B14)  yields 


(B14) 


c  V  5x  .  •  5y  .  ,  Scu  ,  ,  _  Sf2 

S  =  X^  +  Y'^+h^ - l-h  cos  I  ^ 

OS  OS  OS  OS 


(B15) 
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Let  J  =  ^  ^  ^  substitutions  from  table  I  equations 

(1-144),  (1-151),  (1-158),  and  (1-163)  may  be  introduced  to  obtain 

-  =  ^  (l  -  e^  cos^E)+  a(sin  E)(l  -  e  cos  E)fe  ^  +  a 

g  os  \  Os  os/  ('32,6) 

where  E  is  the  eccentric  anomaly. 

The  time  derivative  of  Kepler's  equation  (eqs.  (1-47)  and  (1-52)) 
is 

E(l  -  e  cos  E)  =  n  =  pl/2a“(2/2)  (B17) 

v/here  n  is  the  mean  angular  motion  in  the  orbit.  Eliminate  E  in 
equation  (B16)  using  equation  (B17)  to  obtain 

J  =  pi/2jll/2(i  +  e  cos  E)  II  +  sin  E  ||  +  aV2  sin  E  |^ 

which  can  be  written  as 


J  = 


^  aV2(i  _  e  cos  E)  ||  -  |  a"^^/“^(E  -  e  sin  E) 


5a 


+  i  al/h3 


+  e  cos 


E) 


|E  ,,  1  ^-(1/2) 

OS  4 


(3E  + 


sin 


E) 


oa 


+  ab2 


sin  E 


J  =  - 


,1/2 


2ia 


A 

3s 


[a3/2(E 


e  sin 


E)’ 


.1/2 


3 


fal/' 


(3E  + 


e  sin 


E)J 


ui/2  3g  ^ 

2a  3^  3^ 


(B18) 


where  H  =  i  pl/2al/2(3E  +  e  sin  E)  and 

G  =  aV2(E  -  e  sin  E)  =  n3?/‘^{t  -  t^)  =  pl/2(t  -  t^)  .  Thus, 

3g  i/2 


(B19) 
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Introduce  equations  (B18)  and  (B19)  into  equation  (B15)  to  obtain 


^  ,  SH  ,  ,  acD  ,  ,  ^  Sa 

b  =  -  - - - 1-  +  h  +  h  cos  I 

2a  as  os  os  os 


(B20a) 


A  similar  derivation  will  show  that  Q  is  given  by 


U  ^  ,  SH  ,  ,  Sod  ,  ,  ^  Sfi 

- - - h  ^  +  h  v f-  h  cos  I  V- 

2a  oq  oq  oq  oq 


(B20b) 


Using  equations  (B20)^  equation  (B3)  for  Lagrange's  brackets  may  now  be 
vjritten  as 

r„^l_S/  sS  /  p\  S/q\S,  x  ScDSh  ShSoD 

^  Sit  ^  V  ^  ^Si  -  57  Si 


Sil  S(h  cos  I)  5(h  cos  l)  ^ 

57  5q  Ss  oq 


(B21) 


Using  Jacobian  notation,  the  general  expression  for  a  Lagrangian  bracket 
as  obtained  by  rewriting  equation  (B2l)  is 


-t 


C_  1  \  P’  2a/  (a),h)  (a,h  cos  I) 

-  (s,J  +  — tyi) — 


(B22) 


Substitution  of  h  =  Vpa(l-e^  into  equation  (B22)  gives  the  expres¬ 
sion  from  which  the  Lagrangian  brackets  will  be  evaluated  for  the  chosen 
elements 


I 


f 
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APPENDIX  C 


DISTUEBING  FUNCTIONS  OF  THE  ELEMENTS  IN  TERMS  OF 

COMPONENTS  OF  THE  DISTURBING  ACCELERATION 

By  the  application  of  spherical  trigonometry  it  is  seen  that  the 
components  A^,  A^  of  the  disturhing  acceleration  referred  to  OX^ 

OY ,  and  OZ  in  figure  1  are  given  in  terms  of  C,  and  R  hy 

=  R(cc,j^  cos  V  +  CL2  sin  v)  +  C(a2  cos  v  -  sin  v)  +  W  sin  H  sin  I 

(Cla) 

Ay  =  R(|33_  cos  V  +  pg  -  Pp  sin  v)  -  W  cos  Q  sin  I 

(Clh) 

A^  =  RCVp  cos  V  +  Fg  sin  v)  +  C(T2  cos  v  -  sin  v)  +  W  cos  I 

(Clc) 


where 


ap  =  cos  CO  cos  -  sin  co  sin  ^  cos  I 
ag  =  -sin  CO  cos  Q.  -  cos  co  sin  Q  cos  I 
(3p  =  cos  CO  sin  SI  +  sin  co  cos  Q  cos  I 
Pg  =  -sin  CO  sin  +  cos  co  cos  Q,  cos  I 
Tp  =  sin  CO  sin  I 
Tg  =  cos  CO  sin  I 


(C2) 


Geometrically,  ap ,  Pp,  yp  are  the  direction  cosines  of  the  pericenter 
radius,  and  ag ^  Bg ’  Yz  direction  cosines  of  the  circumferential 

direction  at  pericenter  referred  to  OX,  OY,  OZ,  respectively. 

To  obtain  Ox/ds,  x  -»  y,z  where  s  is  any  of  the  orbital  elements, 
the  coordinates  of  the  object  must  be  expressed  in  terms  of  the  orbital 
elements.  The  (x,y,z)  coordinates  in  the  system  OX,  OY,  OZ  may  be  ex¬ 
pressed  in  terms  of  the  orbital  elements  and  true  anomaly  v  as 
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X 


a(l  -  e2)  . 

1  +  e  cos  V  ^  +  <^2  sin  v 


(C3a) 


y 


a(l  -  e^) 

1  +  e  cos  V 


(PjL  cos  V  +  Pg  sin  v) 


z 


a(l  -  e^) 

1  +  e  cos  V 


(Ti  cos  V  +  Tg  sin  v) 


(031) 

(C3c) 


when  equations  (1-140),  (1-147),  (1-154),  (l-72),  and  (1-104)  from  table 
I  are  used. 


Determination  of  Dg^,  Dg,  D.j. 


From  equations  (C2)  it  is  seen  that  the  derivatives  of  a,  p,  and 
T  with  respect  to  the  elements  a,  e,  and  tp  are  zero.  Equations 

(C3)  are  not  explicit  in  terms  of  the  chosen  set  of  orbital  elements, 
for  sin  v  and  cos  v  are  seen  to  be  functions  of  a,  e,  and  tp  by 

the  following  form  of  Kepler's  equation  obtained  from  equations  (1-18), 
(1-24),  (1-45),  (1-47),  and  (l-52)  of  table  I: 


(t  -  tp)^.l/VV2  = 


tan' 


sin  V 


e  +  cos  V 


e 


-y/l  -  e^  sin  V 
1  +  e  cos  V 


(C4) 


Because  of  the  transcendental  form  of  equation  (C4),  explicit  expressions 
for  sin  v  and  cos  v  in  terms  of  a,  e,  and  tp  cannot  be  obtained. 

However,  the  required  partial  derivatives  may  be  obtained  according  to 
the  following  relations  if  v  is  regarded  as  an  auxiliary  variable: 


3x  rSx'l  Bx  3v 

^  -  y.z  (C5) 


for^  s  equal  to  a,  e,  tp  where  [Bx/ds]  in  brackets  Indicates  that  the 
derivative  is  now  taken  only  so  far  as  s  appears  explicitly  in  equa¬ 
tions  (C3) .  The  derivatives  from  equations  (C3)  are 


Bx 


a(l  -  e^ 


sin  V  -  a2(e  + 


cos  V 


>]- 


x,a  ^  y,P;  z,r 
(C6) 


(l  +  e  cos  v) 

rBx*l  1  -  e^  , 

“  1  +  e  cos  V  ^“l  ^  °^2  (C7) 
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3x  _  a[2e  +  (cos  v)(l  +  e^)] 
(l  +  e  cos  v)2 


(a-(_  cos  V  +  ag  sin  v)^  x,a  ^  y^P;  z,r 

(C8) 


and  from  equation  (C4)  the  derivatives  are 


(C9) 


Bv  ^  _  3  [ 

3a  2'Y^f;/^  ?\3 


(l  +  e  cos  v)^(t  -  t  ) 

5(1  -  e2)-^  ^ 


3v  sin  V 


oe 


1  -  e'^ 


(2  + 


e  cos  V 


§i_  -  rz 

^  “  Va3( 


(l  +  e  cos  v) ‘ 


Hi  -  e2)2 

Combining  equations  (C6)  to  (C12)  according  to  equation  (C5)  yields 
Sx  1  -  e2  ^ 

_  - 2Qg  y  +  „  g3_j^ 

da  1  +  e  cos  v  -l  2 


(CIO) 

(Cll) 

(C12) 


x,a  -  y,P;  z,r  (C13a) 


dx 
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a.{_Ze  +  (cos  v)(l  +  e2)j 
(l  +  e  cos  v)^ 


(a^  cos  V  +  ag  sin  v) 


-  a  sin  V 


2  +  e  cos  V 
(l  +  e  cos  v)2 


ai 


sin  V 


agCe  +  cos  v)j  , 


x^a  ^  y^P;  z,r  (C13b) 


dt^ 


Va(l  -  e2)  [' 


|a^  sin  V  -  UgCe  +  cos  v)  ,  x^a  y^P;  z,T 


(C13c) 


2021- 
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The  desired  results  for  Dg_.  Dg,  and  are  obtained  by  combining 

equations  (Cl)  with  equations  (C13a)^  (C13b)^  and  (C13c)^  respectively, 
according  to  equation  (I3a) .  Rote  that  the  following  conrnon  factors  in 
the  process  have  the  values  shown: 

Coefficients  of  tenns  containing  R: 

(a^  cos  V  tttg  sin  v)^  +  (|3^  cos  v  +  pg  sin  v)^  +  cos  v+ Tg  sin  v)^  =  1 

(a^  cos  V  +  ttg  sin  ’v)  |aq  sin  v  -  a^Ce  +  cos  v)J 

+  (P-,  cos  V  +  Pg  sin  ■v)^Pq  sin  v  - 

+  (l^  cos  V  +  Tg  sin  “''')  sin  v  -  T^ie  +  cos  v)j  = 


-e  sin  V 


Coefficients  of  terms  containing  C: 

(ag  cos  V  -  aq  sin  v) (aq  cos  v  +  ag  sin  v) 

+  (Pg  cos  V  -  Pq  sin  v)(pq  cos  v  +  pg  sin  v) 

+  (Vg  cos  v  -  Tq  sin  v) (fq  cos  v  +  Tp  sin  v)  =  0 
(cl2  cos  V  -  aq  sin  ■v)jjxq  sin  v  -  agCe  +  cos  v)J 


+  (pg  cos  V  -  Pq  sin  v) 
+  (Tg  cos  V  -  Tq  sin  v) 


P-l  sin  V  -  Pp(e  +  cos  v' 


Tq  sin  V  -  TgCe  cos  v, 


=  -(l  +  e  cos  v) 


Coefficients  of  terms  containing  W: 

(aq  cos  V  +  ag  sin  v)sin  S2  sin  I  -  cos  SI  sin  I 


•  (pq  cos  V  +  Pg  sin  v)  +  (y  cos  v  +  sin  v)cos  1=0 


Jaq  sin  V  -  a2(e  +  cos  v)J  sin  sin  I  -  cos  Q.  sin  I 

■  P 


Sin  V 


P2(e  +  cos  v)J  +  |Tq  sin  v  -  1=0 
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The  simplified  results  then  are  obtained  as 


= 


|_x  +  e  cos  V  2  "W  a3(i  _  e^) 


v; 


(t  -  t.  )e  sin  V 


R 


a'^(l  -  6“' 


(t  -  t  )(l  +  e  cos  v)  C  (C14) 


Dg  ==  ( -a  cos  V 


)R  +  a(sin  v)^- 


2  +  e  cos  ~v 
+  e  cos  V, 


(CIS) 


D, 


“\a(l  -  e^) 


—  [(e  sin  v)R  +  (l  +  e  cos  v)c]  (C16) 


Determination  of  Dj^^  Dj 

The  disturbing  functions  of  cd,  K,  and  I  are  obtained  by  combining 
the  derivatives  of  equations  (C3)  with  equations  (Cl)  according  to  equa¬ 
tion  (l3a)  after  equations  (C2)  are  used  to  eliminate  a,  3^  and  T  in 
equations  (Cl)  and  (C3) .  As  this  work  is  rather  long  but  very  direct,  it 
is  omitted.  Results  are: 


=  Cr  (C17) 

Dj  =  Wr  sin  u  (ClS) 

D^  =  Cr  cos  I  -  V/r  cos  u  sin  I  (C19) 
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TABLE  I.  -  EQUATIONS  FOR  THE  TWO -BODY  ORBIT 
(a)  Classification  of  orUits^ 


Type  of 
orbit 

Energy, 

Eg 

Angiilar 

momentum, 

h 

Eccentricity^ 

e 

Semimajor 

axis, 

a 

Semilatus 

rectum, 

P 

Circle 

Eg  =  -u/2r 

h  =  rV 

e  =  0 

a  =  r 

p  =  r 

Ellipse 

<  0 

/  0 

<  1 

>  0 

>  0 

Degenerate 

ellipse 

<  0 

=  0 

=  1 

>  0 

=  0 

Parabola 

=  0 

/  0 

=  1 

=  CO 

>  0 

Degenerate 

parabola 

=  0 

=  0 

=  1 

=  CO 

=  0 

Hyperbola 

>  0 

i  0 

>  1 

<  0 

>  0 

Degenerate 

hyperbola 

>  0 

=  0 

=  1 

<  0 

=  0 

^The  degenerate  conics  included  for  completeness  are  straight 
lines. 


The  twD-hodj^-  equations  describing  motion  on  a  conic  section  may  be 
written  in  myriad  forms.  The  following  particularly  useful  equations 
have  been  selected  for  tabiuLation  without  derivation. 

It  is  assumed  thau  e  a.nd  r  are  never  negative.  Special  cases 
and  restrictions  necessary  for  real  equations  are  indicated.  It  is 
also  obvious  that  many  of  the  equations  become  indeterminate  or  un¬ 
bounded  for  certain  conditions.  (Symbols  are  defined  in  appendix  A.) 


TABLE  I.  -  Continued.  EQUATICaiS  FOR  THE  TWO-BODY  ORBIT 
(b)  Continued.  Two-body  eq.uatlons 


0 


1-42: 


1-43: 


Orbital  inclination 


tan  I 


0  <  I  <  n 


Bin  I 


1-44:  cos  I  «3  ^ 

h 


Mean  anomaly 


1-45: 

M  =  n(t  -  tp). 

e  <  1 

1-46: 

“  v(t  -  tp). 

e  >  1 

1-47: 

=  E  -  e  sin  E, 

e  <  1 

1-48: 

a  -F  +  e  sinh 

e  >  1 

1-49: 

M  =  n 

Mean  angular  motion 

1-50 1 

M 

e  <  1 

X  -p 

1-51: 

- 

2« 

P  ' 

e  <  1 

1-52: 

e  <  1 

1-53: 

-  e^)^,  e  <  1 

1-54: 

V  = 

-  1)'^,  e>l 

"5  * 

1-55: 

-in 

Semilatus  rectum 


1-56:  p  B  a(l  -  e^) 

1-57:  “  2*(l  +  e  cos  v) 

I-58I  = 

1-59:  = 

1-60:  w  (l  +  +  2e  cos  v) 

yZ 


1-61: 

1-62: 


=  ^(1  - 

1  -  e  cos  E^ 

2(ii/r)  -  V2 


e  <  1 


Semilatus  rectum  (concluded^ 

1-641 

P 

Zi-pTa 
”  Tp  +  ra 

1-65: 

=  r(l  +  cos  v),  e  =  1 

1-66: 

s=  2rp,  e  =»  1 

1-67: 

=  r,  e  =  0 

Period 

1-68: 

P  =  — ,  e  <  1 

n 

1-69: 

=  2jt^2,  e  <  1 

Radius 

1-70: 

P 

1  +  e  COB  V 

1-71: 

■y/x®  +  y2  +  z2 

1-72: 

a{l  -  e^) 

1  +  e  cos  V 

1-73: 

- 

a(l  -  e  cos  E)^  e  <  1 

1-74: 

= 

a(l  -  fe  cosh  F)^  e  >  1 

1-75: 

+  tan2  1^,  e  1 

I-76i 

= 

P/  e  «  0 

1-77: 

r  = 

7  sin  ^ 

1-78: 

xx  +  yy  +  zz 

r 

1-79: 

= 

“  *  I  =  0  =  n 

r 

1-80: 

= 

t/—  e  Bin  V 

fP 

1-81: 

r  B 

ue  cos  V 

rS 

Radius  at  pericenter 

1-82: 

= _ E- 

P  1  +  e 

1-83: 

=  Ea  -  ra 

1-84: 

=  a(l  -  e) 

1-85: 

=  |,  e  =  1 

Radius  at  apocenter 

1-63  J 


=  4^  e  <  1 


1-86: 

1-87: 


^e. 


P 

1  -  e 
2a  -  Tp 


I- 


a(l  +  e) 


1-140: 

1-141: 

1-142: 

X  =  r  cos  V 

1-143: 

w  x(cos  m  cos  £2 

1-144: 

=  a(cos  E  -  e). 

1-145: 

=  a(cosh  F  -  e)> 

1-146: 

=  rp(l  -  tan^  |) 

1-147: 

1-148: 

1-149: 

Y  =  r  sin  V 

1-150: 

=  -x(sin  (0  cos  fl  • 

1-151; 

=  a  sin  E, 

1-152; 

=  -aVe^  -  1  sinh 

1-153: 

=  2rp  tan  < 

1-154; 

TABLE  I.  -  Continued.  EQUATIOITS  FOR  THE  TWO-BODY  ORBIT 
(L)  Continued.  Two-body  equations 

Rectangular  coordinates 
X  =  r(cos  u  cos  SI  -  sin  u  sin  SI  cos  I) 

=  r  cos  u,  I  =  0  = 


y  =  r(cos  u  sin  SI  +  sin  u  cos  SI  cos  I) 
=  r  sin  u,  I  =  0  =  Q 


x(sln  (u  006  n  +  coa  m  sin  a  cos  I)  -  y(8ln  to  sin  a  -  cos  a>  cos  a  cos  I)  +  y  cos  m  sin  I 


z  »  r  sin  u  sin  I 
Velocity  components 


^  ^  cos  I(e  coa  a>  +  cos  u)  +  cos  Sl(e  sin  m  +  sin  u)] 

=  "  ^  (e  sin  (jj  +  sin  u),  I  s=  0  =* 

X  =  -  sin  V 


=  -aE  sin  E,  e  <  1 

=  aF  sinh  F,  e  >  1 


y  =  [cos  SI  cos  I(e  cos  o)  +  cos  u)  -  sin  Sl{e  sin  o  +  sin  u)] 

=  (e  cos  a>  +  cos  u)^  I  «  0  =  12 

Y  «  (e  +  COB  v) 

=  a"^!  -  e®  E  cos  E,  e  <  1 
=  -a'^e^  -  1  p  cosh  e  >  1 

z  =  sin  l(e  cos  co  +  cos  u) 


TABLE  I.  -  Concluded.  EQUATION’S  FOR  THE  TWD-BODY  ORBIT 


(■b)  Concluded.  Two-body  equations 


Path  angle 

1-166: 

tan  i)r  =  — 

rv 

1-167: 

xi  +  yy  +  zz 
h 

1-168: 

e  sin  V 
~  1  +  e  cos  V 

1-169: 

er 

=  —  sin  V 

P 

1-170: 

sin  i|f  = 

V 

1-171: 

XX  +  yy  +  zz 
rV 

1-172: 

_  e  sin  V 

Vl  +  e^  +  2e  cos  v 

I-]73: 

cos  =  £Z 

V 

1-174: 

h 

“  rV 

1-175: 

_  1  +  e  cos  V 

\/l  +  e^  +  2e  cos  v 

Argument  of  perlcenter 

1-176: 

CQ  U  -  V 

Ascending  node 

1-177: 

tan  n  =  -^r— 

1-178: 

sin  a  =  - X - 

1-179: 

r. 

cos  Q  - 

■l/h^  + 

T  X  y 
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TABLE  II.  -  DEHIVAICnVES  OF  ORBITAL  ELEMEHTS  AED  PARAMETERS®' 

(a)  Time  derivatives  of  orbital  elements  due  to 
perturbations 


E-1202 
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TABLE  XX.  -  Continued.  DERIVATIVES  OF  ORBXTAL  ELEMENTS  AND  PARAMETERS 


(a)  Continued.  Time  derivatives  of  orbital  elements  due  to 
perturbations 


Time  of  perl center 


-  [» ^  -  f)  -  -  *p)  ?]=}' 

■  vVJ  {”*)  -  =<*  -  *  (^'•} 


0  ^  e^l 


0  5^  e  5^  1 


-  2e  cos  V  +  £2£J!;  _  3  +  -  +  3e  ml.  0  <  e  <  1 


sin  V  -  3  -2  _®=-_  T 

M  vr^ 


=  £  JP  (—I  fe  +  3e  +  +  cos  v)\ 

V  ^  ^  |l  _  e2  \e  1  +  e  cos  V  / 


+  cos  vjN|,  0  <  e  <  1 

The  results  for  tp  for  e  >  1  are  Identical  to  the  preceding  results  for 
0  <  e  <  1,  but  vlth  E  replaced  by  IF  and  Vl  -  e2  replaced  by  1  "  !• 

2  r 

t  =  - - - ^  I  r2(sln  v)(cos2v  +  3  cos  v  +  l)]c 

5|a(l  +  cos  v)^  '• 

-  [2  cos^v  +  4  cos2v  +  cos  v  -  e  =  1 

=  P  §[[(4  cos2v  +  7  cos  v  +  .4)sln  v]  ^  ^  ^  ^  1 

v(l  +  cos  v)  ■y  n  I  ^  5(1  +  cos  v)  J  J 


1  r  /  6t2  +  m2 


iN  +  3vT 


2  yS  I  4qi2  ^  jjp 


1 ,  p  - 


2  ^  y2  LcS  +  r2 


w  +  r/RC  -  CR 

^  vl  2  2 

''Uc'^  +  r'^ 
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TABLE  II.  -  Concluded.  DERIVATIVES  OF  ORBITAL  ELEMENTS  AND  PARAMETERS 

(b)  Time  derivatives  of  orbital  parameters  due  to  both 
orbital  motion  and  perturbations 
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TABLE  III.  -  COMPOLSMTS  OF  THE  DISTURBING  ACCELERATION 


E-1202 
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I 


4Y 


TABLE  lY.  -  QUALITATIVE  EFFECTS  OP  THE  DISTURBIIiG  ACCELERATIOM  COMPONENTS® 


rr 


Derivative 

of 

Coir.j^cnent 

parameter 

T 

tl 

C 

R 

W 

Semlmajor 
axlSy  a 

Always  + 

0 

Always  + 

1’ 

and  2’ ,  + 

0 

3* 

and  4 ’ ,  - 

Seir.llatus 
rectbuny  p 

Always  + 

1*  and  2’j 
3’  and  4’ , 

,  + 

Always  + 

0 

0 

Radius  of 
perloenter 

’  '■p 

Always  + 

1'  and  2' , 
3'  and  4’, 

,  + 

Always  + 

V 

3» 

and  2‘ ,  - 
and  4S  + 

0 

Eccentrlcl ty , 

,  e 

j  1  and  4,  + 

1'  and  2' , 

- 

Vl  -  e^  -  1 

COS  V  ■>  -i-i - = - i.  +  1 

1^ 

and  2 ’ ,  + 

1 

2  and  3,  - 

3'  and  4* , 

+ 

cos  V  <  ~  - 1— - 

1 

3' 

and  4* ,  - 

0 

Mean  anomaly, 

,  M 

- ’  and  2' ,  - 

:  1  ’  and  4  ’ , 

- 

'  1 ’  and  c ‘ ,  - 

GCfi 

+ 

3’  and  4’,  + 

2’  and  3* , 

+ 

3'  and  4’ ,  + 

Vl  +  beS  -  1 

'  ''  <  2e 

“ 

Tr;e  anomaly. 

V 

1’  and  2’ ,  - 

1"  and  4", 

_ 

1 ’  and  2 • ,  - 

1* 

3*  and  4’,  + 

2”  and  3", 

+ 

3 ’  and  4 ’ ,  + 

2* 

and  3’ ^  - 

0 

Argument  of 

1»  and  2’ ,  + 

1"  and  4", 

and  4' ,  - 
and  3’ ,  - 

I  and  II,  - 
III  and  IV, 

perlcenter. 

OJ 

3’  and  4',  + 

2"  and  3", 

+ 

3’  and  4’ ,  + 

2» 

Ascending 

I  and  II,  + 

node,  ii 

0 

0 

III  and  IV, 

Inclination, 

I 

0 

• 

I  and  lY,  + 

0 

II  and  III, 

It  muat  be  ren.erabered  that  the  table  applies  only  to  perturbation  ter"J  of  the  clerlvSt^^r  RMu^^ta  Sje 
opposite  In  every  case  for  negative  disturbing  components.  derivative.  Results  are 


gure  1.  -  DiagraTTi  showing  various  orbital  parameters  if.  generei; 
OX;  OY;  OZ  system  and  location  of  0X2_,  OXg;  and  OX3  coordinate 

axes  in  alternate  systems. 


W»»»wrr: 
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Figure  1.  -  Diagram  showing  various  orbital  parameters  in  , 
OX,  OY,  OZ  system  and  location  of  OX^^,  OXg,  and  OX3  coor- 
axes  in  alternate  systems . 
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